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Abstract 

The higher-order perturbative corrections, beyond leading logarithmic accuracy, to the BFKL evo¬ 
lution in QCD at high energy are well known to suffer from a severe lack-of-convergence problem, 
due to radiative corrections enhanced by double collinear logarithms. Via an explicit calculation of 
Feynman graphs in light cone (time-ordered) perturbation theory, we show that the corrections en¬ 
hanced by double logarithms (either energy-collinear, or double collinear) are associated with soft 
gluon emissions which are strictly ordered in lifetime. These corrections can be resummed to all 
orders by solving an evolution equation which is non-local in rapidity. This equation can be equiv¬ 
alently rewritten in local form, but with modified kernel and initial conditions, which resum double 
collinear logs to all orders. We extend this resummation to the next-to-leading order BFKL and 
BK equations. The first numerical studies of the collinearly-improved BK equation demonstrate 
the essential role of the resummation in both stabilizing and slowing down the evolution. 

Keywords: QCD, Renormalization Group, Color Glass Condensate, Hadronic Collisions 
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1. Introduction 

It is by now well established that the Balitsky-JIMWLK hierarch}!^ [IH2] and its mean field ap¬ 
proximation known as the Balitsky-Kovchegov (BK) equation |3] govern the high-energy evolution 
of scattering amplitudes in presence of non-linear effects (multiple scattering and gluon saturation) 
responsible for unitarization. Some of the most remarkable recent developments in that context re¬ 
fer to the first calculations of the next-to-leading order (NLO) corrections [5H7] to the B-JIMWLK 
and BK equations. These new developments parallel and extend previous efforts, towards the end 
of nineties, which established the NLO version mm of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation [MS] — the linearized version of the BK equation which applies so long as the scat¬ 
tering is weak. Although the BFKL and B-JIMWLK equations are based on a common evolution 
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mechanism, they differ in the way how they treat the scattering problem: the BFKL equation deals 
only with single scattering, as appropriate for a dilute target, whereas the B-JIMWLK hierarchy 
includes the interplay between evolution and multiple scatterings. The former is usually written 
in transverse momentum space, as an equation for the unintegrated gluon distribution, while the 
latter is formulated in terms of transverse coordinates (better suited suited for implementing the 
eikonal approximation) and keeps trace of the multiple scattering of the individual partons in the 
projectile — each of them represented by a Wilson line. Such differences explain the difficulty 
to adapt to the NLO B-JIMWLK evolution the ‘collinear resummations’ originally developed in 
the context of NLO BFKL |13H16| . which aim at improving the convergence of the perturbative 
expansion for the BFKL kernel. 

The collinear resummations refer to perturbative corrections, starting at NLO, which are en¬ 
hanced by large, single or double, transverse logarithms. Without a proper resummation, which, 
strictly speaking, goes beyond the order-by-order expansion of the BFKL kernel, these large loga¬ 
rithms deprive the NLO BFKL formalism of its predictive power. 

There is no reason to expect this lack-of-convergence problem to be attenuated by the non¬ 
linear terms in the B-JIMWLK equations; indeed, the ‘collinear’ corrections arise from regions 
in phase-space where the scattering is weak and the non-linear effects are negligible. This was 
anticipated in a semi-analytic study m and later on confirmed by the numerical observation that 
adding a unitarity constraint (in the form of a ‘saturation boundary’) to the NLO BFKL equation 
does not help improving the stability of the solution [18] . Very recently, while our work was being 
completed, this has been corroborated by a numerical study [19] of the NLO BK equation [5| : the 
numerical solution turns out to be unstable (the scattering amplitude decreases with increasing 
energy and can even turn negative) for the physically interesting initial conditions. As also shown 
in Ref. |19j . this instability can be traced back to a large NLO correction to the BFKL kernel 
enhanced by a double transverse logarithm. This kind of correction, which can be associated with 
the choice of the reference scale in the energy logarithm, is well understood at BFKL level, where 
it is successfully resummed to all orders by the schemes proposed in Refs. |13H16j . It is our main 
objective in this paper to propose a similar resummation at the level of the BK equation. 

More precisely, our goals are twofold: first, we would like to unambiguously identify the origin 
of the double-collinear logarithms in Feynman graphs to all orders and devise a method for their 
resummation; second, we would like to reformulate this resummation as a change in the kernel of 
the BK equation, which is energy independent. Unlike the corresponding method in the context of 
NLO BFKL, where the resummation is generally implemented in double Mellin spac{0 [T3] ll4|. our 
resummation will be directly implemented in transverse coordinate space, in order to be consistent 
with the non-linear structure of the BK equation. 

Concerning the first objective above, our main finding is that the double-collinear logs arise 
due to a reduction in the longitudinal phase-space for the high-energy evolution, as introduced by 
the condition that successive gluon emissions be strictly ordered in lifetime. The interplay between 
this ‘kinematical constraint’ and the double transverse logarithms has already been recognized 
in the literature m [20] (see |21) for a recent discussion and more references), but we are not 
aware of any systematic derivation of this prescription from Feynman graphs. To emphasize that 
this is indeed non-trivial, we notice that double collinear logs are also generated by diagrams 
with anti-time ordering, but they mutually cancel when all such graphs are summed together (see 


^Note however some similarity between our strategy and that proposed in [16], where the cu-shift in Mellin space 
has been approximately reformulated as an improvement of the BFKL kernel in transverse momentum space. 
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the discussion in Sect. below). This observation helps understanding the peculiar way how the 
double transverse logs arise in the context of the NLO BK calculation in [5]. The main outcome 
of this diagrammatic analysis is Eq. which governs the evolution in the double-logarithmic 
approximation (DLA): it resums to all orders the perturbative corrections in which each power of 
the coupling is accompanied by a double logarithm (either energy-collinear, or double collinear). 

Eq. ( |17[ ) however is non-local in ‘rapidity’ (the logarithm of the longitudinal momentum, which 
is our evolution variable), so it does not fully match our goals for a collinearly-improved evolution 
equatiorj^ To cope with that, in Sect. we demonstrate that the non-local equation 0 can be 
reformulated in a local form, modulo an analytic continuation and a reshuffling of the perturbative 
expansion. The new, local, equation ( |30[ ) involves an ‘improved’ kernel and (for consistency) a 
modihed initial condition, which both resum double-collinear logs to all orders. 

It is then straightforward to extend this resummation to the BEKL and BK equations and 
thus obtain the collinearly-improved BK equation (32), which is our main result in this paper. 
It is furthermore possible to promote this result to full NLL accuracy, by adding the remaining 
NLO BK corrections from Ref. [^. Notice however that the NLO terms include single transverse 
logarithms, which may require additional resummations, as was already the case in the context of 
NLO BEKL [IMS]. 

Finally, in Sect. we present the first numerical studies of the resummed BK equation (32). 
These studies clearly demonstrate the role of the resummation in both stabilizing and significantly 
slowing down the evolution: the saturation exponent extracted from the numerical solution is 
smaller by, roughly, a factor of two than in the absence of the resummation. 


2. The double-logarithmic limit of the BFKL equation 


In order to fix the notations and for comparison with the more refined results that we shall 
later obtain, it is instructive to recall the derivation of the ‘naive DLA’, by which we mean the 
version of this approximation which neglects the time-ordering of successive gluon emissions, from 
the leading-order (LO) BFKL equation |10m2j . The LO BFKL equation resums the perturbative 
corrections in which each power of the QCD coupling as = agNc/n, assumed to be fixed and small, 
is accompanied by the energy logarithm Y = \t\.{s/Qq) (the ‘rapidity’), with s the center-of-mass 
energy squared and Qq a characteristic transverse momentum scale introduced by the target. In 
this high-energy leading-log approximation (LLA), valid when oigY > 1, it is consistent to treat the 
scattering and the evolution in the eikonal approximation. The LO BFKL equation can then be 
written as the linearized version of the BK equation mm, f.e. as an equation for the high-energy 
evolution of the scattering amplitude TxyiY) of a quark-antiquark dipole, with a quark leg at 
transverse coordinate x and an antiquark leg at transverse coordinate y, which undergoes weak 
scattering off a generic target (a nucleus, or a ‘shockwave’): 


dTxyiY) _ a 


dY 


= ^ / d?zMxy.[Ta,z{Y) + T,y{Y)-TxyiY)] . 


This equation involves the ‘dipole’ version of the BFKL kernel, 
M = __ 


( 1 ) 

( 2 ) 


® Collinearly-improved versions of the BK equation which are non-local in rapidity have been proposed too in the 
literature |21[ I22| . but they suffer from some shortcomings, concerning either the systematics of the resummation 
(for the approach in [22]), or its feasibility in practice (for |21p. 
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Figure 1: Diagrams illustrating one step in the BFKL evolution of a qq dipole via the emission of a soft ‘real’ gluon 
(the gluon fluctuation lives at the time of scattering). The target is represented as a shockwave. In the (eight) 
corresponding ‘virtual’ graphs, the gluon line is not crossing the shockwave. 


which describes the emission of a soft gluon with transverse coordinate 2 ; by either the quark or 
the antiquark leg of the dipole, followed by its reabsorption (see Fig. [^. In the limit of a large 
number of colors Nc —>■ 00 , the positive quantity {a.s/27r)M.xyzd'^z can be interpreted [23] as the 
differential probability for the splitting of the original color dipole {x, y) into a pair of dipoles 
[x, z) and (z, y). The first two terms within the square brackets, Txz and T^y, are the ‘real’ terms 
describing the scattering of the daughter dipoles, whereas the last one, Txy, is the ‘virtual’ term 
expressing the reduction in the probability for the parent dipole to survive at the time of scattering. 

Eq. Q is valid so long as the scattering is weak, T <C 1, meaning that all the dipoles look small 
on the target resolution scale: (a; — regime, the integration over z in the 

r.h.s. of Eq. 0 becomes logarithmic when the daughter dipoles are much larger than the original 
one, i.e. for |a; — z| ~ \z — y\ r = \x — y\. Eor such configurations, the fast decreases of the 
dipole kernel, Mxyz — Iirt, — z^^ , is partially compensated by the rapid increase of the scattering 
amplitudes for the daughter dipoles: Txz — Tzy oc {x — z)^. Note that this property, namely that 
the dipole scattering amplitude grows, roughly, like the area of the dipole in the transverse plane, 
is indeed satisfied by the initial condition at low energy E ~ 0 and is also preserved by the high 
energy evolution in the (double-logarithmic) approximation discussed belovQ 

The same arguments also imply that, in the region where the z~integration is logarithmic, the 
‘virtual’ term Txy oc is much smaller than the ‘real’ ones and can be neglected. Its only effect 
will be to cut off the logarithmic phase-space at small dipole sizes \x — z\ ~ r, or \y — z\ ~ r. 
Indeed, when z ^ x or z ^ y, the linear combination in the r.h.s. of Eq. Q vanishes due to 
‘real’-vs.-‘virtual’ cancellations and to ‘color transparency’ {T(Y,r) —)• 0 as r —)• 0). 

To write down the double-logarithmic version of the BEKL equation, it is convenient to factor 
out the strong r^-dependence of the amplitude and write Txy{Y) = r'^Q'^AxyiX). Eocusing on 
the logarithmically-enhanced contributions only, we can then average over azimuthal angles and 
impact parameters, in order to replace AxyiY) —)• A{Y,r'^). The ensuing equation for A{Y,r‘^), in 
the (naive) double logarithmic accuracy, is most conveniently written in integral form, as 


A{Y, r^) = ^(0, r^) as 



"i/Qo dz^ 


^A(Yuz^ 


(3) 


where the upper limit l/Qg is truly the unitarity limit, i.e. the transverse scale where the scattering 
becomes strong: T{Y,r) ~ 1 when r ~ l/Qo- This integral equation can be solved via iterations: 
A = Yl'^=o ^ with of order d”. Eor instance, for the simple initial condition ^(0,r^) = 1, 


■^In the full BK equation, for small enough r^, the same argument applies to a large phase-space 1/Qs S> (x — z)^ 
with Qs the saturation momentum in the target. 
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one finds (with p = In l/(r^(5o)) 


A^\Y,p) = asYp, A^^\y,p) 


a^Y'^p^ 

4 


•••, a^^\y,p) 


{asYpY 

( n !)2 


(4) 


whose sum is recognized as a modified Bessel function, A{Y,p) = lQ{2YasYp). 

As obvious from the previous discussion, the ‘naive DLA’ equation Q resums terms of the 
type {agY p)” to all orders. This is a common limit of the BFKL and DGLAP evolution, which 
however is not very useful in practice, since valid only in a very narrow regime (namely, when 
OigYp > 1, but such that oisY <C 1 and oisP 'C 1). In the next section, we shall devise a more 
general double-logarithmic approximation, which resums all the perturbative corrections where ds 
is enhanced by exactly two logarithms — either Yp oi p^. 


3. From Feynman graphs to the DLA equation 

To obtain a more general DLA formalism, we shall perform a diagrammatic calculation, using 
light-cone (time-ordered) perturbation theory, of two successive steps in the high energy evolution 
of the dipole scattering amplitude. That is, we shall consider 2-gluon graphs in which one of the 
emitted gluons is much softer (in the sense of having a much smaller longitudinal momentum k~^) 
than the other one. We use conventions where the dipole projectile is an energetic right mover 
and we work in the light-cone gauge A~^ = 0. The use of ‘old-fashioned’ time-ordered perturbation 
theory is particularly convenient for our purposes, in that it allows for an economical and physically 
transparent classification and evaluation of the relevant Feynman graphs in the approximations of 
interest. 

A rather compact way to organize the calculation has been devised in Ref. |24j : successive gluon 
emissions by the projectile, which are time-ordered and also strongly ordered in k~^, are generated 
by repeatedly applying an evolution ‘Hamiltonian’ on the 5-matrix for the high-energy scattering 
between the projectile and the target color field A~. This 5-matrix is built from Wilson lines (one 
for each parton in the projectile), whose number keeps increasing in the course of the evolution, 
due to additional gluon emissions. For the original dipole, one has Txy = 1 — Sxy with 

Sxy = ^ tT^[V^Vy] , where Fj = Pexp ji^ y dx+A"(a;+, ®)t“| . (5) 

The Hamiltonian acts via functional differentiation w.r.t. A~ and describes the emission of a 
soft gluon out of any of the preexisting Wilson lines, followed by its reabsorption (by either the 
same or a different Wilson line). After all the emissions have been produced by acting with the 
Hamiltonian, one must integrate over all the emission times and average over the background field 
A~ to construct the physical scattering amplitude. This last step though (the functional average 
over A~) is irrelevant for the present purposes and will be left unspecified in what follows. 

The first two steps in this evolution generate 2-gluon graphs like those shown in Fig. where 
the gluon with longitudinal momentum p'^ is emitted first and it is harder than the gluon k'^ 
(p"*" > k'^). The topologies (in the sense of time-orderings) selected in Fig. are quite special, 
in that they contribute already to llaQ : in Fig. [^a, the hard gluon (p"*") is real and is emitted 

®The time-orderings contributing to LLA can be divided in two classes: (i) if the hard gluon is ‘real’ (it crosses the 
shockwave), then it is emitted before the soft one (ti < ti), and reabsorbed after it (t 2 > T 2 )\ (ii) if the hard gluon 
is ‘virtual’, then the 2 gluons have no overlap in time with each other — the hard gluon is emitted and reabsorbed 
either before the emission of the soft gluon (fi <t 2 < vi) or after the absorption of the latter {t 2 > t\ > T 2 )- 
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Figure 2: Diagrams with two gluons which are ordered in longitudinal momentum (p~^ > k'^) and also in lifetime 
(Tp > Tk); (a) a real-real graph; (b) a virtual-real graph. 


before the soft one (k~^), but reabsorbed after it; in Fig. @b, the hard gluon is virtual and it is 
both emitted and reabsorbed prior to the emission of the soft gluon, which is real. Beyond LLA, 
other time orderings become important as well and will be later considered (see Fig. [^. 

We shall first evaluate the 2-real-gluon graph in Fig. [^a. After integrating over all emission 
times, within the ranges — oo < < ri < 0 and 0 < T 2 < t 2 < oo, one finds the following 

contribution to the change in dipole S'-matro® (below, = f d^w and f^ = f 


f O O q 


^ipiu-x) ^ip-(x-u) ^ik-(z-y) ^ik-{u-z) 


■(z-v)Jk-(u-z) P'P 


ppkk 

X 


dA;+ n* 


dp"* 


P 


p-f 


'4 




lk+ W p+ + k+l2 p+ + k+ 


(p-fc)^ 

k^ 


( 6 ) 


In the integrals over k'^ and p"*", the upper limit is the longitudinal momentum of the quark 
and antiquark in the original dipole, while the lower limit is the longitudinal scale at which 
the scattering probes the dipole wavefunction; that is, the overall rapidity interval available to the 
evolution of the projectile wavefunction is y = ln(q'*“/g(|). The denominators in the second line 
come from time integrations and can be recognized as the usual ‘energy’ (here, in the sense of p~) 
denominators of light-cone perturbation theory. For instance, 


-b A :+P + ^ Tp-\- Tk 

where Tp = 2p'^ jp^ = l/p~ is the lifetime of the hard gluon fluctuation, as determined by the un¬ 
certainty principle, and similarly for r^. The integral over p'^ is logarithmic provided p'^ dominates 
both energy denominators, that is, so long a^p"'' > k~^{p^/k‘^), or Tp > Tk- Hence, to leading 
logarithmic accuracy for the longitudinal logarithm, one can replace Tpjijp -|- Tk) — &{Tp — Tk)- 
In the BFKL regime, one assumes that there is no strong hierarchy between the transverse 
momenta, |fc| ~ |p|, so the condition Tp > Tk is automatically satisfied when p'^ > k^- In that 


®To keep expressions simple, we use the large-W limit at intermediate steps, but some of the final results, notably 
the DLA equation (171, are valid for any Nc- 

^For the purposes of power counting, one can use |fc| ~ |fe| and |p| ~ |/5 — fe|; indeed, the difference between e.g. 
k and k is due to the scattering off the target, which is a comparatively small effect in the high transverse momenta 
(or small dipole sizes) regime of interest. 
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case, one can freely integrate over transverse momenta in expressions like Eq. (©, to generate the 
Weizsacker-Williams propagators of the soft gluons, according to 


/ 



P P ip-{x—z) 

\ O O ^ 


27r {x — 2)2 


i a:* — 2 ® 


( 8 ) 


After also summing over all possible connections for the two emitted gluons, one builds the relevant 


product of dipole kernels (i.e., M.xyu-M.uyz for the sequence of emissions illustrated in Fig. [^. 

However, this is strictly correct only so long as the transverse phase-space is by itself not 
logarithmic, meaning so long as E ^ p, where p = measures the logarithmic separation 

in transverse scales between the original dipole, with size r = 1/Q, and the target, with size 
l/Qo- In the end, the transverse integrations in Eq. Q are restricted to this range, e.g. Qq < 
(see below). For sufficiently large values of p, one opens the phase-space for a logarithmic 
integration over p^, which favors relatively large values |p| \k\. In this regime, the theta-function 

Q{Tp — Tk) = &{p'^ — k'^{p^/k‘^)) becomes relevant and its effect is to reduce the longitudinal phase- 
space, roughly from Y toY — p. 

To the accuracy of interest, i.e. to correctly keep both the corrections of orders agY and 
generated when integrating out the hard gluon p '^, the constraint Tp > can be enforced directly 
in coordinate space, like p'^u'^ > k'^z^. Here, we have anticipated that the corrections of the form 
otsP^ come from emissions which are strongly ordered in transverse sizes, such that the daughter 
dipoles are much larger than the parent one. In this regime, 

|2 — ~ |2 — y| ~ |2 — » \u — x\ zz \u — y\ ^ r = \x — y\ , (9) 

and u refers to any of the sizes, \u — x\ or \u — y\, of the first pair of daughter dipoles, while 2 
similarly refers to the daughter dipoles produced by the second splitting After performing the 
momentum integrals in Eq. Q, summing over all the possible connections for both emitted gluons, 
and adding the other splitting sequence (where the gluon at 2 is emitted from the dipole {x,u)), 
one finds the following result from the 32 time-ordered graphs with two ‘real’ gluons (at large Nc). 



Except for the theta-function enforcing time-ordering, this is recognized as the effect of two con¬ 
secutive steps in the LO BFKL evolution. 


To this result, one must add contributions coming from virtual graphs, evaluated to the same 
accuracy. The ‘real-virtual’ graphs in which the harder gluon (p"*") is virtual, whereas the softer 
one {k~^) is real, are the only ones that matter to the accuracy of interest. Consider first the 32 
such graphs whose topologies (i.e. time-orderings) exist already at LLA, namely those where the 
two gluons have no overlap in time with each other (an example is shown in Fig. [^b). They give 



( 11 ) 


®When integrating over generic values u and z, like in Eq. (101 below, one can set u = max {\u — x\, \u — y\) and 


2 = max(|z - a:|, \z - y\, \z - u\). 
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In the BFKL context, this contribution is used to regulate the short-distance singularities of Eq. ( |10[ ) 
asu ^ X and it —)■ y at a scale set by the original dipole size: u> r. In the present context, it plays 


a similar role (as anticipated in Eq. Q), except for the fact that only the time-ordered piece of © 
is needed for that purpose. That is, albeit the virtual graphs included in Eq. 0 do not naturally 
involve any time ordering, it is nevertheless useful to distinguish between the respective time- 
ordered (TO) and anti-time-ordered (ATO) contributions, by inserting 1 = ©(xp — r*,) -|- ©(r^ — Tp) 
in the integrand of Eq. ©• (Here and from now on, Tp = p'^u^ and = k^z"^.) Then the 


TO piece must be combined with the 2-real-gluon contribution in Eq. (10), which is itself time- 
ordered, whereas the ATO piece is to be considered together with other virtual-real graphs, which 
are naturally ATO and will be discussed below. 

From now on, we shall limit ourselves to the strict double-logarithmic approximation (DLA), 
where each power of cxs is accompanied by either Yp or p^. The corresponding contribution of 
Eq. (10) can be isolated by taking the single scattering approximation and restricting the integra¬ 
tions over u and 2 : according to Eq. (§. This allows for simplifications like 


■^^uyz-^^xyu — 




1 _ Q Q Q ~ T 

>^xu^uz^zy — 


+ T,y ~ 2Tiz^) 


( 12 ) 


For subsequent discussions, it is important to stress that, to DLA, it is only the last emitted gluon 
(the one with the largest tranvese size z) which contributes to scattering. Then the integrals over 
p'^ and u are both logarithmic, as anticipated, and can be evaluated as 



(13) 


where the logarithmic variables Y = /k'^) and p = ln( 2 ^/r^) refer to the phase-space available 

to the hard gluon Note that we have implicitly assumed above that Y > p, so that the integral 
over p~^ has indeed support for any u > r. This can be recognized as the condition for the lifetime 
Tk = k~^z'^ of the soft gluon fluctuation be (much) smaller than the ‘lifetime’ Tq = of the 

original dipole (the duration of the quantum process which has produced that dipole, e.g. the 
fluctuation of the virtual photon in DIS). 

To summarize, by integrating out the intermediate gluon p~^, one has produced, besides the 
expected LLA contribution OgYp, also a contribution which can be interpreted as a NL© 

correction to the BFKL kernel for the emission of the soft gluon k~^. This correction matches 
the respective piece (that enhanced by a double transverse logarithm) of the full NL© result in 
Ref. [3]. The last remark might suggest that the remaining 2-gluon graphs, that have not been 
considered so far and which correspond to other time orderings, do not contribute to order OsP^- 
But this is not quite true: contributions of this order arise from all the diagrams which are anti¬ 
time-ordered (AT©), in the sense that the lifetime of the hard gluon is shorter than that of the soft 
one (to DLA, at least). Topologically, the class includes two types of diagrams: (i) real-virtual 
graphs where the hard gluon is virtual and overlaps in time with the soft gluon which is real (some 
examples are the graphs la, lb, 2a, 3a, 3b, 4a, and 4b in Fig. [^; (ii) real-real graphs where the 
hard gluon is emitted after, and absorbed before, the soft one (see graph 2b in Fig. |^. To this 
genuinely ATO diagrams, one must add the AT© pieces of the virtual-real graphs without overlap 
in time (see graphs Ic, Id, 3c, and 3d in Fig. which represent the ATO part of graphs like that 
in Fig. 0b. left over from the earlier calculations), to cancel UV divergences and introduce an 
effective short-distance cutoff equal to r (cf. the discussion after Eq. 0). 
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Figure 3: Pattern of cancellations (to DLA) in 2-gluon graphs with anti-time-ordering (rp < Tk). 




When evaluating graphs of the type (i) and (ii) above mentioned, one finds that the time 
integrations over the overlapping region produce a factor like 


P 

p~ + k~ 


Tk 

Tp + Tk 


~ 0(rfe 


Tp), 


(14) 


where the theta-function approximation in the r.h.s. holds in the double-logarithmic region. This 
theta-function cuts off the rapidity phase-space at the scale p (with p <Y) and thus produces a 
contribution oc as anticipated: 



It turns out however that all the terms of order asP^ generated by these ATO graphs exactly 
cancel each other. These cancellations can be understood as either the cancellation of ‘infrared’ 
(large z^) logarithms between self-energy and vertex corrections, or as real vs. virtual cancellations 
for hard gluons whose scattering is not measured at DLA (cf. the discussion after Eq. (dll))- For 
instance, the combinations of graphs la and lb in Fig. or 3a and 3b, belong to the first category, 
whereas 2a and 2b belong to the second. The sum of la and lb leaves an uncompensated UV 
divergence, which is regulated at the scale r after also adding the ATO pieces of Ic and Id (and 
similarly for 3a, 3b, 3c and 3d). Finally, graphs 4a and 4b mutually cancel because the color current 
associated with the hard gluon (and responsible for the emission of the soft gluon) has a different 
sign in 4a as compared to 4b. Interestingly, the order-dsp^ contribution of the real-real graph 2b is 
such that it would exactly compensate the respective contribution of all the time-ordered real-real 
graphs previously discussed. This explains why, when the calculation is organized in such a way 
that the real-real graphs are all grouped together, as in Ref. [5] (which used the standard Feynman 
rules in momentum space), the net effect of order cksP^ is rather seen to arise from the sum of the 
virtual-real graphs. 

The above cancellation pattern for the ATO graphs naturally generalizes to higher orders, i.e. 
to graphs involving an arbitrary number of strongly-ordered (in p~^) gluon emissions. This leads us 
to the main conclusion in this section, namely the fact that the net contributions to DLA come fully 


9 














































































from graphs which, within time-ordered perturbation theory, have exactly the same topology as the 
graphs contributing to LLA, but with the additional constraint that the successive gluon emissions 
must be strictly ordered in lifetimes. This implies that the perturbative corrections enhanced by 
the double logarithms Yp or can be resummed to all orders by solving a modified version of the 
DLA equation Q, which includes the time-ordering condition: 


A{q^,r^) = A{0,r‘^) + as 


ri/Ql ^ rq- 

Jr^ Jq+ 


0 + 


dfc+ 

fc+ 


A(A:+,, 


(16) 


In what follows, we shall mostly use logarithmic variables, with the target scales and Qq being 
the reference scales: Y = ln(g+/( 7 (]'), Yi = ln(A;+/g(|), p = ln(l/r^(5Q), pi = \n{l/ z^Qq). Then 
Eq. (16) becomes 


A{Y, p) = A(0, p) + as / dpi 


rY-p+pi 


dTiA(yi,pi). 


(17) 


where a step function Q(Y — p + pi) is implicitly assumed within the integrand, to ensure that 
Ti > 0. This integral equation determines the function A(Y, p) for all positive values of Y and p, 
but the most interesting physical regime lies at T > p. 


4. Resummed kernel for DLA, BFKL, and BK evolutions 

As compared to its ‘naive’ version in Eq. (§, the DLA equation with time ordering 0 is 
non-local in rapidity, as it can be best appreciated by rewriting it as a differential equation for the 
T-evolution. This non-locality complicates the practical applications and, more importantly for 
our present purposes, it makes it quite tricky to extend this equation beyond DLA accuracy (albeit 
similar non-local versions of the BEKL and BK equations have been proposed in the literature; see 
Ref. [21] for a recent discussion). In what follows, we shall construct an alternative version of this 
equation, which is equivalent to 0 in the most interesting physical regime at K > p and which 
is local in rapidity. The generalization of this local equation to BEKL and BK will then be almost 
straightforward, up to NLL accuracy. 

At the mathematical level, it is more convenient to first solve the following problem: 


fiY,p) = f{0,p)+as dpi 


rY-p+pi 


dYif(Yi,pi), with /(0,p) = (5(p). 


(18) 


Given its solution f{Y,p), one can immediately construct the solution to Eq. 0 for arbitrary 
initial conditions according to 


A{Y,p)= dpi/(y,p - pi)A(0,pi). 
Jo 


(19) 


It is straightforward to solve Eq. (18) via iterations, to find 

oo _ k k _ 1 

/(y,p) = 5 (p) + ^/(")(y,p), with /W(y,p) = 0 (y-p)^^L^^Z^ 


( 20 ) 


10 






where the 0-function arises from the longitudinal phase-space integration. For the purpose of the 
physical interpretation, it is useful to keep in mind that f{Y, p) is essentially the unintegrated gluon 
distribution in the dipole. The presence of the 0-function in the solution reflects the fact that, 
in order to emit a soft gluon, its lifetime must be smaller than the coherence time 

Tg = IQ'^ of the original dipole. Summing the above series we arrive at the explicit form 


fiY,p) = 6{p) + eiY-p), 


'as{Y - p) 


Ii {2^as{Y- 


( 21 ) 


where Ii is the modified Bessel function. Neglecting for the moment the 0-function one can show 
that the function above admits an integral representation in the complex plane; namely, one has 
/(y, p) = Q{Y — p)f{Y, p), where the new function f{Y, p) is defined as 


t+icxD 




L ^exp 




a.c 


1-C 


(Y-p) + il-Op 


( 22 ) 


for any positive Y and p. This can be viewed as the analytic continuation of the original function 
/(y, p) outside the physical range Y > p. Making the change of variables 7 = ^-|-q;s/( 1 — 0; 
can recast Eq. (22) in the form of a standard Mellin representation, that is 


fiY,p)= i ^•>^(7)exp[a^XDLA(7)d" +(1-7)p] • 


y-ioo 


Here, the “characteristic function” is determined by 
d5XDLA(7) = ^ -(1 - 7) -h V(1 - 7)2 4 as 


CXa 


a? 


+ 


2a3 


( 1 - 7 ) ( 1 - 7 )^ ( 1 - 7 ) 


+ ••• 


(23) 


(24) 


where all the poles at 7 = 1 visible in the r.h.s. are merely an artifact of expanding Xdla( 7) in 
series of as- The resummed answer is clearly hnite at 7 = 1. The Jacobian J{j) induced by the 
change of variables is related to the characteristic function and reads 


J(7) = 1 - a^XDLA(7) = 1 - 


a.s 


+ 


3a 


A,2 


+ 


(25) 


(1_7)2 

The existence of a Mellin representation together with the exponentiation in Y (as manifest 
in the integrand in Eq. (23)) demonstrate that the function f{Y,p) obeys an evolution equation 
which is local in Y. Namely, Eq. (23) is tantamount to the following integral equation 


f{y,p) = f{0,p)+as[ dYi [ dpiJCoi^Aip-Pi)fiYi,pi), 
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IQ 


with the kernel /Cdla(/ 3) defined as the inverse Mellin transform of Xdla( 7), that i: 

dl{ 2 y/Osp"^) asP^ {dsP^f 2n3x 

/Cdla(/5) = — =1-^ H-77,-h 0 {{asp ) ), 




asp2 
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(26) 


(27) 


®Such a kernel (albeit written in momentum space) has already appeared in a work focusing on the collinear 
improvement of the NLO BFKL equation m- In that context, this kernel was obtained as an approximation to the 
equation for the BFKL eigenvalue with cu-shift [Si¬ 


ll 

























(Ji is the Bessel function) and the initial condition f{0,p) obtained as the limit of Eq. (23) at the 
unphysical point Y = 0 : 


/(0,p) = 6{p) - 


(28) 


To summarize, the solution to Eq. (26) with the kernel (27) and the initial condition (28) exists 


for any positive values Y and p. For Y > p \t reduces, by construction, to the original function 


/(y, p) in Eq. (21). The importance of this construction is that it can be immediately generalized 
to the evolution of the dipole amplitude, which can be thus rewritten as a local equation in Y. 
First, we define the analytic continuation of A{Y,p) according to (cf. Eq. ([T^) 


dpi/(y,p - pi)^(0,pi). 

Jo 


(29) 


This new function coincides with the physical amplitude A{Y, p) for Y > p. For general, positive, 


values of Y and p, it obeys an equation similar to Eq. (26), that is, 

A{Y,p) = A{^,p) + as [ dYi [ dpi/C dla(p - Pi)^(yi, Pi), 

Jo Jo 


(30) 


with an initial condition ^4(0,p) which follows from Eqs. (29) and (28). For illustration, consider 
two interesting initial conditions, namely ^(0, p) = 1, which has the advantage of simplicity, 
and ^(0, p) = p, which is the limit of the McLerran-Venugopalan (MV) model for dipole-nucleus 
scattering in the single scattering approximation [25]. One easily hnds 


M(0,p) = < 


)[1 + Jo0] 


P 

2 


1 + Jo (p) + !" Ho (p) Ji (p) - ^ Hi (p) Jo (p) 


for M(0,p) = l, 
for M(0, p) = p. 


(31) 


where we have temporarily used the notation p = 2\/a^iA and where is the Struve function. 


Eq. (30) is the sought-after local version of the DLA equation for the dipole amplitude: for 


Y > p, its solution coincides, by construction, with the respective physical amplitude, i.e. with the 
solution to the non-local equation ©• Notice that this rewriting of the DLA evolution in local 
form is tantamount to a complete reshuffling of the perturbation series: both the kernel in Eq. (30) 
and the initial condition in Eq. resum double-collinear terms of the type (asp^)"" for any n. 
For instance, the very first iteration of this equation generates all the terms linear in dsV, i.e. 
the terms of the type OgYp{^sP^)"‘ with n > 0, that would be produced by iterating the original 
equation ( |17| ) to all orders. Remarkably, even though both the kernel and the initial condition 
exhibit oscillations as functions of p, their combined effect within equations like (30) or ( |26| ) yields 
a solution which is positive definite in the physical region Y > p, order by order in 6is (e.g., this 


produces the perturbative solution (20) for f{Y,p)). 

As we now explain, it is rather straightforward to promote this local DLA equation into a more 
complete equation, which includes the right BFKL and BK physics to NLL accuracy. To that aim, 


and starting with Eq. (30), we shall make backwards the steps leading from the LO BFKL equation 


([^ to the ‘naive’ DLA equation ([^, that is: 

(i) we use the full expression for the dipole scattering amplitude, and more precisely its analytic 
continuation T{Y,p) = q~pA{Y,p) (which coincides with the physical amplitude for Y > p); 
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(ii) we return to the use of transverse coordinates in our notations, meaning that we replace 
p = \n{l/r^Ql), p-pi= ln(zVr2), f{Y, p) = f^y{Y), and 2f{Y, z^) ^ f^,{Y) + f,y{Y)- 

(iii) we restore the full dipole kernel by replacing {r‘^/z^)dz‘^ —)> M-xyz^^lT^'-, 

(iv) we reintroduce the virtual term and at the same time remove the infrared and ultraviolet 
cutoffs on the integral over z, since they are not needed anymore; 

(v) we replace the argument of the (additional) kernel /Cdlaj that is, p — p\ = ln(z^/r^), 
according to ln( 2 :^/r^) —)• \jLxzr^‘yzr-, wher ^^ ^ixzr = ln[(a: — zf'j{x — y)"^] and = {x — y)^. 

The last step above is the only one which is new as compared to the original discussion in 
Sect, [^and will be thoroughly justified in a moment. We are thus led to the following equation 


dfxyjY) 

dY 


a 


2^ / h zAAxyz I^DLA 


;/C, 


LxzrLyzr 


fx.{Y)+fzy{Y)-fxy{Y)-fxz{Y)fzy{Y) , (32) 


where we have also added the non-linear term familiar from the BK equation, to account for 
multiple scattering and thus ensure unitarization. Eq. (32) improves over the LO BK equation by 
resumming the double-collinear logs, i.e. the perturbative corrections of the form (d^p^)’^, to all 
orders. Importantly, this resummation affects both the kernel and the initial conditions. 

So far, the initial condition at K = 0 has been specified only in the weak scattering regime 
where T <C 1 and the precise normalization of T(y = 0,r) was unessential (since the respective 
evolution was linear). For the purposes of the non-linear equation (32), however, we need to 
hx this normalization. To that aim, it is convenient to use the MV model [25], which amounts 
to exponentiating the amplitude for a single scattering: T(0,r) = 1 — where To(^) = 

e~^A{0,p), with p = ln(l/r^(5Q) and M(0, p) as shown in the second line in Eq. (31). 

To fully motivate Eq. (32), we still need to explain our choice for the argument of /Cdla in this 
equation. Clearly, the replacement made at point (v) above is irrelevant at the present interest 
(improved LLA), but has two important virtues: (a) it switches off the effect of the collinear 
resummation for small daughter dipole sizes, \x — z\ < r, or |y — z| < r, where this resummation 
should indeed play no role, and (b) when expanded to second order, i.e. /Cdla(/ 5) 1 — dsp^/2, it 

precisely matches the double-logarithmic term contained in the full NLO BK result |5|. The last 
feature makes it straightforward to formally extend Eq. (32) to full NLL accuracy: it is sufficient 
to add to its r.h.s. all the NLO BK corrections computed in [^, except for the double-log term 
that has already been included in the kernel. 

Let us finally remind that, strictly speaking, the solution to Eq. (32) can be trusted only for 
sufficiently small values of p < K (which in turn requires Y to be large enough, oigY > 1, in order 
to significantly evolve away from the ‘unphysical’ initial condition). In practice though, we expect 
the BFKL evolution encoded in Eq. (32) to eventually wash out the oscillations introduced by the 
initial condition at large p and thus progressively built a physical tail including at p > K. This 
will be checked via numerical calculations in the next section. 


5. Numerical tests 

In this section we present a brief selection of first numerical studies which illustrate some subtle 
issues previously discussed (like the interplay between local and non-local evolution equations) and 
also some physical consequences of the resummation. 


^°Some caution is required here, since there is a small region of integration where the product appearing under 


the square root can be negative; in this case, it is enough to let Lsaz.rLy 2 ,r 
with < 0). 


and Ji —>■ Ii (cf. Eq. (271 
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A(Y,p) 




Figure 4: Left: the DLA solutions for A(Y, p), cf. Eq. ( |17| l (full lines), and respectively A{Y, p), cf. Eq. ( |30| (dashed 
lines), for various rapidities and (physical) initial condition A(0, p) = 1. Right: the characteristic function asXil) 


numerically extracted from the kernel in Eq. (321 vs. its LO (BFKL) and ‘NLO’ approximations. 


Consider first the double-logarithmic approximation. In the left-hand plot of Fig. we show 
the function A(Y,p), which we recall is related to the dipole amplitude T{Y,p), as obtained from 
two different approachej^ (a) by directly solving the non-local equation with initial condition 
A{f),p) = 1 , and (b) by solving the local equation (30) with the initial condition shown in the first 
line of Eq. ( |31[ ) (which, strictly speaking, truly yields the analytic continuation A{Y,p)). The 
respective solutions are supposed to coincide only at p <Y, where they both represent the actual 
physical result. This is indeed confirmed by the numerical simulations. On top of that, for p > Y, 
we see that the (physical) solution to Eq. (0 is independent of p and equal to AiY, p) = cosh y/og Y 
whereas its analytic continuation A(Y, p) shows (non-physical) oscillations which are inherited 
from the initial condition. 


We now move to the collinearly-improved BK equation (|32|). Numerically, we solve the evolution 

By acting 


equation following a strategy similar to the one described in the Appendix of Ref. 
with the kernel in this equation on the power-like test function one can numerically extract 
the characteristic function q;sx( 7 ) (the would-be Mellin transform of the resummed kernel). In 
the right-hand plot of Fig. we compare the function otsxil) thus obtained for the particular 
value ots = 0.25 (black triangles) with the respective LO (BFKL) approximation asXo( 7 ) (red 
squares) and with a ‘NLO’ approximatiorj^ dsXNLo( 7 )) obtained by keeping only the ds term 
in the expansion of /Cola, that is, /Cola —^ K,^no{p) = 1 — OisP^j^. The solid lines show the 
expected analytic results for the LO and NLO curves and their agreement with the numerical 
results is a powerful check that the numerical procedure is under control. As manifest on this 
figure, the behavior near 7 = 1 is strongly influenced by the higher order corrections. This can also 


be understood by inspection of the DLA approximation Xdla( 7 ) in Eq. (24). At NLO accuracy. 


both cases, we found that a simple discretisation of the integral and a Euler method to solve the rapidity 
evolution was sufficient to reach good numerical accuracy. 

^^fn this section, by ‘NLO’ we refer to the inclusion of the NLO corrections which are enhanced by double collinear 
logarithms, finite NLO corrections being neglected. 
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Figure 5: Numerical solutions to the BK equation for the dipole amplitude at strict LO (l.e. Eq. (321 with /Cdla f), 
NLO (meaning with kernel /Cdla —^ A^nlo ), and after resummatlon (l.e. with the full kernel /Cdla of Eq. \27\ ). The 
dashed line In Gg. (c) Indicate the transition between Y < p and Y > p; dotted lines are the direct result of the 
numerical simulation, while solid lines have been matched to the expected physical behaviour for p > Y, l.e. T oc e“^. 


Xdla(7) exhibits a cubic pole at 7 = 1, the second term in the r.h.s. of Eq. (24), with a negative 
residue which makes the function Xnlo( 7 ) unstable in the collinear limit 7 —?■ 1 (in particular, 
there is no saddle point on the real axis). By contrast, the all-order resummation ensures a smooth 


behavior near 7 = 1, as already noticed after Eq. (24). For as = 0.25, the function 7 ( 7 ) is seen to 
be almost flat for 7 > 0.5. 


A crude estimate of the saturation lin^^ based on the DLA result in Eq. (21) yields 127! 


Vn 


with A, = 


da. 


1 + Aas 


(33) 


which is significantly smaller than the respective LO result (no resummation) Abfkl — d.SSds |25) . 
This suggests that the reduction of the longitudinal phase-space coming from time-ordering and 
giving rise to collinear double logs leads to a considerable reduction in the speed of the evolution. 

This expectation is indeed confirmed by the numerical solutions to Eq. (32). In Fig we show 
the results for = 0-25 and for an initial condition of the MV type, with Al(0, p) = 1 (and hence 
Al(0, p) as given in the first line of Eq. (m). As before, the results with all-order resummation (cf. 
Fig [ 5 ] c) are compared to the respective predictions of LO BFKL (cf. Fig[^a) and to the ‘NLO’ 
results obtained by using /Cnlo(p) = 1 ~ ds/9^/2 (cf. Figj^b). The latter are highly unstable and 
physically meaningless — the evolution rapidly leads to a negative scattering amplitude — as it 
could have been anticipated in view of the pathological behavior of the corresponding characteristic 
function Xnlo( 7) “ Fig- 0 Similar instabilities have been recently observed m in numerical 
simulations of the full NLO BK equation and they have been traced back to the large double- 
logarithmic terms ~ ctsP^ in the NLO kernel, in agreement with our present findings. By contrast, 
the evolution with the fully resummed kernel, shown in Fig[^c, is perfectly smooth. We also see 
in Fig l^c that the non-physical oscillations at p > Y introduced by resummation in the initial 
condition tend to disappear at larger rapidities. Finally, by comparing the LO results in Fig[^a to 
the resummed ones in Fig[^c, one clearly sees the anticipated reduction in the evolution speed. A 
quick estimate of the saturation exponent from the numerical results in Fig[^c yields A^ ~ 0.55, 
in remarkable agreement with the crude DLA estimate in Eq. (33). 


^We recall the saturation line Ps{Y) is defined by the condition that T{Y,p) ~ 1 when p — ps{Y). 
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